Anomalous and spin Hall effects in hep cobalt from GGA + U calculations 
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We have calculated the intrinsic anomalous and spin Hall conductivities, spin and orbital mag- 
netic moments and also spin polarization of Hall currents in hexagonal cobalt within the density 
functional theory with the generalized gradient approximation (GGA) plus on-site Coulomb inter- 
action (GGA+Z7). The accurate all-electron full-potential linearized augmented plane wave method 
is used. We find that the anomalous Hall conductivity (AHC) (cr^y) and orbital magnetic moment 
(m ) of cobalt with the magnetization being along the c-axis (m//c) calculated in the GGA+U 
scheme with U = 1.6 eV and J = 0.9 eV agree rather well with the corresponding experimental 
values while, in contrast, the a^ y and m from the GGA calculations are significantly smaller than 
the measured ones. This suggests that moderate U = 1.6 eV and J = 0.9 eV are appropriate for Co 
metals. The calculated AHC and spin Hall conductivity (SHC) (o~ S ) are highly anisotropic, and the 
ratio of the Hall conductivity for m//c to that for the magnetization being in-plane (m//ab) from 
the GGA+C/ calculations, is ~16.0 for the AHC and ~6.0 for the SHC. For m//c, the spin-up and 
spin-down Hall currents are found to flow in the same direction. The magnitude of the calculated 
Hall current spin polarization (P H ) is large and P H — —0.61. Remarkably, for m//ab, the spin-up 
and spin-down Hall currents are predicted to flow in the opposite directions. This indicates that 
the Hall current contains both spin-polarized charge current and pure spin current, resulting in the 
magnitude of the Hall current spin polarization (P H — —1.59) being larger than 1.0. 

PACS numbers: 71. 15. Mb, 71.70.Ej, 72.15.Gd, 72.25.Ba 



INTRODUCTION 



Spin-dependent electronic transports have attracted 
intensive attention recently because they not only are in- 
teresting from the view-point of fundamental physics but 
also have fascinating technological applications. Anoma- 
lous Hall effect (AHE), discovered in 1881 by Hall 1 , is an 
archetypal spin-related transport phenomena and hence 
has received renewed interests in recent yearsi 2 - The ordi- 
nary Hall effect in nonmagnetic conductors is driven by 
the Lorentz force due to the applied magnetic field while 
the AHE depends on magnetization in ferromagnetic 
materials 3 . The dipole magnetic field produced by the 
magnetization is too small to explain the observed AHE 
in ferromagnets, and thus it has been realized for long 
time that the AHE must be caused by the spin-orbit in- 
teraction (SOI). Therefore, several competing SOI-based 
mechanisms have been proposed. Extrinsic mechanisms 
of skew scattering^ and side jump£ refer to the modi- 
fied impurity scattering induced by the SOI, and hence 
should in principle depend on the nature and density of 
the impurities. Surprisingly, extrinsic side jump contri- 
bution to the anomalous Hall conductivity (AHC) turns 
out to be fairly independent of either scattering strength 
or defect densityi 2 ^ Another mechanism arises from the 
anomalous velocity of the Bloch electrons caused by the 
SOI, discovered by Karplus and Luttinger— , and is thus 
of intrinsic nature. Interestingly, this intrinsic mecha- 
nism has recently been reinterpreted in terms of the Berry 
curvature of the occupied Bloch states^ - — Furthermore, 
recent first-principles studies based on the Berry phase 



formalism showed that the intrinsic AHE is important 
in various materialsi 10 ' 11 In particular, in itinerant fer- 
romagnets such as Fe, the intrinsic AHC given by first- 
principles density functional calculations with the gener- 
alized gradient approximation (GGA)^ has been found 
to agree rather well with the experimental AH C 12 ' 13 . 

Nonetheless, the physical origin of the AHE in cobalt 
and some other ferromagnets are still not fully under- 
stood. Recent GGA calculations predicted an intrinsic 
AHC of ~480 S/cm in hep Co with the magnetization 
along the c-axis, being smaller than the corresponding 
experimental value by about 40 %>M Density functional 
calculations with either the local density approximation 
(LDA) or GGA have been rather successful in describ- 
ing many physical properties such as crystal structure, 
elastic constant and spin magnetic moment, of itinerant 
ferromagnets Fe, Co and Ni (see, e.g., Ref. [IH and ref- 
erences therein). However, LDA and GGA calculations 
often fail in describing relativistic SOI-induced phenom- 
ena in these itinerant magnets. For example, the the- 
oretical values of orbital magnetic moment account for 
only about 50 % of the measured ones in Fe and Co 16 i 17 
and the calculated magnetocrystalline anisotropy energy 
of hep Co is even wrong in sigr>i&. This failure of the 
LDA/GGA is generally attributed to their incomplete 
treatment of the 3d electron-electron correlation in these 
systems. Several theoretical methods that go beyond the 
LDA and GGA, such as the orbital-polarization correc- 
tion (OPC)^ and LDA/GGA plus on-site Coulomb inter- 
action U (LDA/GGA+J/)^"— schemes, have been devel- 
oped for better description of the SOI-induced phenom- 
ena in magnetic solids. Indeed, the orbital-polarization 
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correction has been found to bring the calculated or- 
bital moments in itinerant magnets such as Fe and Co 
in good agreement with experiments^ 2 -. It has also been 
demonstrated that the correct easy axes and the mag- 
nitudes of the magnetocrystalline anisotropy energy of 
Co^ 3 - as well as Fe and Ni^i can be obtained within either 
the OPC— or LDA+[7 scheme 2 ^. Furthermore, recent 
GGA+U calculations 25 showed that including appropri- 
ate on-site Coulomb interaction U could reduce the the- 
oretical AHC by half, bringing the calculated intrinsic 
AHC in good agreement with the latest experiments. 

The AHE is closely connected to the spin Hall effect 
(SHE)— ~— , which refers to the transverse spin current 
generation in nonmagnetic materials by an applied elec- 
tric field. In analogous to the AHE, the mechanisms for 
the SHE can be separated into intrinsic and extrinsic 
types. Both intrinsic AHE and SHE are caused by the 
opposite anomalous velocities experienced by the spin- 
up and spin-down electrons when they move through 
the relativistic energy bands in solids under the elec- 
tric field. In nonmagnetic materials where the numbers 
of the spin-up and spin-down electrons are equal, these 
opposite transverse currents would give rise to a pure 
spin current. On the other hand, in ferromagnets, where 
an unbalance of spin-up and spin-down electrons exists, 
the same process would result in a charge current. The 
SHE has attracted intensive current interests both ex- 
perimentally and theoretically since the theoretical pro- 
posal for its intrinsic mechanism in semiconductor s 27 ' 28 , 
because it would enable us to control spins without mag- 
netic field or magnetic materials. Indeed, first-principles 
density functional calculations on the intrinsic spin Hall 
conductivity (SHC) have recently carried out for both 
nonmagnetic semiconductor s 30 ' 31 and metals 3 - 1 - - — . And 
the calculated SHC o 32 ' 33 are found to agree well with the 
corresponding experimental values in many metals such 
as A12 5 ., Pt^&i^ 7 ., Au and Pd— , suggesting that the in- 
trinsic SHE is dominant in the high resistivity regime 
in pure metals. Nonetheless, recent experimenta l 39 ' 40 
and theoretical 4 ^— studies showed that certain impu- 
rities could also give rise to large SHE in impure metals. 
However, the SHE in ferromagnets remains unexplored. 

In this work, we therefore carry out GGA+J7 calcula- 
tions of intrinsic anomalous and spin Hall conductivities 
and relativistic band structure as well as spin and orbital 
magnetic moments of hep Co. The primary objective 
of this work is to better understand the AHE and its 
anisotropy in hep Co. In particular, we would examine 
how the theoretical AHC and related physical quantities 
would be affected when the on-site Coulomb interaction 
is taken into account in the GGA+U scheme. Secondly, 
we would study the intrinsic SHE and its anisotropy in 
hep Co. Among the elemental ferromagnetic transition 
metals, cobalt is promising for future spintronic devices 
because of its large magnetocrystalline anisotropy energy 
and highly anisotropic AHC. A knowledge of the spin- 
polarization of the intrinsic Hall current is important for 
the spintronic applications. Therefore, another objective 



of this work is to evaluate the degree of spin-polarization 
of the charge Hall current by using the calculated AHC 
and SHC. This paper is organized as follows. In the 
next section, we briefly describe how the intrinsic Hall 
conductivities are calculated within the linear-response 
Kubo formalism as well as the numerical method and 
computational details used here. In Sec. Ill, we present 
the calculated intrinsic Hall conductivities, magnetic mo- 
ments, relativistic band structure and also the spin po- 
larization of charge Hall currents. We then compare our 
results with available experiments and also previous cal- 
culations. In Sec. IV, we summarize the main conclu- 
sions drawn from the present work. 

II. THEORY AND COMPUTATIONAL 
METHOD 

An intrinsic Hall conductivity of a solid can be eval- 
uated by using the Kubo formalism^ 6 -. The intrin- 
sic Hall effect comes from the static limit (uj—0) of 
the corresponding off-diagonal element of the optical 
conductivit y 25 ' 30 ' 46 . Following the procedure for the cal- 
culations of the magneto-optical conductivity^ 7 -, we first 
calculate the imaginary part of the off-diagonal elements 
of the optical conductivity 

4 2) M = -JW E E - /knOMWIkn') 

x (kn|^|kn')]5(7iw-e„'J 

(1) 

where V c is the unit cell volume, huj is the photon energy, 
|kn) is the nth Bloch state with crystal momentum k, Vj 
is the j component of the velocity operator, and e n ' n — 
e kn ' — eu_ n . Here ji is the i component of the current 
operator which is — evi and j{<Jk,Vi} for the anomalous 
and spin Hall effects, respectively^ - And also, k ^ i ^ j. 
We then obtain the real part from the imaginary part by 
a Kramers-Kroning transformation 

m 2 f°° ,u/er (2) (w') 

a\co) = -P > (2) 

T Jo OJ Z — 0J Z 

where P denotes the principal value of the integral. The 
intrinsic Hall conductivity is the static limit of the off- 
diagonal element of the optical conductivity a^\uj = 0). 
We notice that the anomalous Hall conductivity of bec 
Fo 10 ' 47 and also the spin Hall conductivity of fee P t 32 ' 48 
calculated in this way are in good agreement with that 
calculated directly by accounting for the Berry phase cor- 
rection to the group velocity. 

Since all the intrinsic Hall effects are caused by the 
spin-orbit interaction, theoretical calculations must be 
based on a relativistic band theory. Here the relativistic 
band structure of hep Co is calculated using the highly 



accurate all-electron full-potential linearized augmented 
plane wave (FLAPW) method, as implemented in the 
WIEN2K code^ 9 . The self-consistent electronic struc- 
ture calculations are based on the GGA for the exchange 
correlation potential^. To further take (i-electron cor- 
relation into account, we introduce on-site Coulomb in- 
teraction U in the GGA+C/ approac h 19 ' 20 . The dou- 
ble counting correction scheme proposed by Czyzyk and 
SawatzkySi is used here. U = 1.6 eV and J = 0.9 eV, 
which were found to give the correct orbital magnetic 
moment for hep Co, are used (see Sec. Ill below). Fur- 
ther calculations with a larger U value of 1.9 eV and 
2.5 eV are also performed to examine how the variation 
of U may affect the calculated anomalous and spin Hall 
conductivities. 

The experimental lattice constants^ 6 - a=2.51 (A) and 
c=4.07 (A) are used here. The muffin-tin sphere radius 
(Rmt) used is 2.2 a.u. The wave function, charge den- 
sity, and potential were expanded in terms of the spher- 
ical harmonics inside the muffin-tin spheres and the cut- 
off angular moment (L max ) used is 10, 6 and 6, respec- 
tively. The wave function outside the muffin-tin sphere 
was expanded in terms of the augmented plane waves 
(APWs) and a large number of APWs (about 80 APWs 
per atom, i. e., the maximum size of the crystal mo- 
mentum K max = 8/ R m t) were included in the present 
calculations. The improved tetrahedron method is used 
for the Brillouin-zone (BZ) integration 51 . To obtain ac- 
curate ground state charge density as well as spin and or- 
bital magnetic moments, a fine 69x69x37 grid of 176157 
fc-points in the first BZ was used. 
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FIG. 1. (color online) Calculated anomalous Hall conductiv- 
ity ofj with and without on-site Coulomb interaction U, for 
the magnetization being either along the c-axis (a) or in the 
a-b plane (b), as a function of the inverse of the number of 
fc-points Nk- 



III. RESULTS AND DISCUSSION 

A. Anomalous and spin Hall conductivities 

Calculated anomalous and spin Hall conductivities 
with and without on-site Coulomb interaction are listed 
in Table HI Like the calculation of the magnetocrys- 
talline anisotropy energy of bulk magnet o 16 ' 17 , a very 
fine fc-point mesh is needed for the anomalous and spin 
Hall conductivity calculations. Therefore, we perform 
the Hall conductivity calculations using several very 
fine fc-point meshes with the finest fc-point mesh being 
137x137x73. The calculated anomalous Hall conductiv- 
ity <7y is plotted as a function of the inverse of the num- 
ber (Nk) of k points in the first Brillouin zone in Fig. 
1. The calculated values of o{\ arc fitted to a first-order 
polynomial to get the converged theoretical (i. e., the 
extrapolated value of of; at Nk—oo) (see Fig. 1). The 
theoretical anomalous Hall conductivity a{j and also spin 
Hall conductivity afj obtained in this manner are listed 
in Table H 

First of all, we notice that taking the on-site Coulomb 
interaction into account via the GGA+f/ scheme affects 
the calculated AHC af, and SHC of, significantly (see 



Table [J and Fig. [TJ. In particular, when the magnetiza- 
tion is along the c-axis (m//c), the theoretical a xy from 
the present GGA calculations is 484 S/cm. When the 
on-site Coulomb interaction is taken into account in the 
GGA+J7 scheme, the calculated AHC a xy increases sig- 
nificantly with increasing U. For example, at U = 1.6 eV, 
the theoretical a xy becomes 618 S/cm, being increased 
by ^30 %. If U is further increased to 2.5 eV, the cal- 
culated a xy is 827 S/cm. Similar behavior is found for 
the a yz with the magnetization being in the a-b plane 
(m//ab) (see Table Q] and Fig. [lj. Secondly, as reported 
before in Ref. [H], the intrinsic AHC in hep Co is highly 
anisotropic. The AHC for m//c is more than 10 times 
larger than that for m/ / ab (see Table HJ . 

The theoretical SHCs of nonmagnetic hep metals were 
recently reported--. Interestingly, we find that the cal- 
culated <rf y for hep Co is larger than that of any 3d hep 
transition metals (see TableUand Ref. l34h . Furthermore, 
hep Zn was found to have the largest anisotropy among 
the nonmagnetic metals studied in Ref. H3- The calcu- 
lated SHC of hep Co is rather anisotropic. The SHC for 
m//c is more than 6 times larger than that for raj jab 
(Table HI, and hence the anisotropy is nearly the same as 
that of hep Zn. The calculated a yz is largest among all 
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TABLE I. Calculated anomalous [afj (S/cm)] and spin [afj 
(ftS/ecm)] Hall conductivities as well as spin [m s (/is/atom)] 
and orbital [m (/xb /atom)] magnetic moments of hep Co 
with magnetization being along either the c-axis or in the a-b 
plane. The on-site Coulomb interaction U — 1.6, 1.9 and 2.5 
eV are used and J = 0.9 eV is kept the same. The estimated 
experimental value of the scattering independent Hall con- 
ductivity alj d as well as previous theoretical intrinsic af and 
extrinsic side jump conductivity afj are also listed for com- 
parison. Note that the estimated al" d should be compared 
with the theoretical af plus theoretical af- 1 . 
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the 3d hep transition metals. 

The AHC af y from the present calculations is in 
good agreement with previous theoretical calculations 14 , 
whilst, in contrast, the AHC a yz presented here is signif- 
icantly smaller than that reported in Ref. (see Table 
|T|). The precise reasons for this discernable discrepancy 
in a yz between the present and previous calculations 1 ^ 
are not known. Note that a yz is one order of magnitude 
smaller than af y (see Fig. 1). On the one hand, the dis- 
crepancy could be attributed to the fact that much more 
fc-points could be included in the BZ integration in Ref. 
Il4l where much more efficient Wannier function interpo- 
lation technique was used. Nevertheless, Fig. 1 indicates 
that the theoretical a yz values presented in Table I arc 
converged with respect to the number of fc-points used. 
In Ref. [ill it was pointed out that the calculated afj 
could be very sensitive to the fine details of the band 
structure especially near the Fermi level. As Figs. 2 and 
3 indicate, the band structure of hep Co is rather com- 
plex near the Fermi level. On the other hand, therefore, 
the discrepancy could also be caused by the fact that the 
band structure of hep Co generated by the Wannier func- 
tions might not precisely reproduce the one calculated 



by using plane wave pseudopotential method in Ref. Il4 . 
Another reason could be the fact that the dipole matrix 
elements for Eq. (1) are calculated using the full APW 
wave functions in the present calculations while they were 
calculated using Wannier functions in Ref. [l4]- 

In Ref. [14], the experimental scattering-independent 
AHC was derived from formula a" ld ~ pf y /p X x using 
the old measured Hall and longitudinal resistivities from 
Refs. l52land l53l respectively, and is 813 S/cm (see Table 
U) . This experimental AHC value should contain both the 
intrinsic AHC af y and extrinsic scattering-independent 
side jump contribution af y . The theoretical side jump 
af. y from the recent ab initio calculations^ is 217 S/cm. 
Therefore, for U = 1.6 eV, the theoretical af plus 
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217 S/cm is almost in perfect agreement with 
the estimated experimental value of a" y d . Nevertheless, 
it should be cautioned that this perfect agreement could 
be fortuitous because in the estimation^ of the experi- 
mental a™ y the Hall and longitudinal resistivities were 
taken from two separated old measurements^^. Fur- 
thermore, the theoretical side jump contributions listed 
in Table U were calculated by considering only the short- 
range disorder in the weak scattering limit^ 4 " and hence 
may not reflect the total side jump mechanism. As will 
be reported in the next subsection, the calculated spin 
and orbital magnetic moments for U = 1.6 eV are also in 
good agreement with the experiments, indicating that the 
appropriate U value for Co metals should be about 1.6 
eV. The experimental AHC a y " d was also estimated using 
more complicated formula a y n z d = af y (p^ z / pf y )(p xx / p zz ) 
in Ref. [lj, and is 150 S/cm. The reported theoretical 
valued of the side jump a yz is -30 S/cm. Consequently, 
the theoretical value (8.9 S/cm at U = 1.6 eV) of a yz 
plus a yz is much too small in comparison with the es- 
timated experimental a % yz . Again, given the fact that 
Pyz/Pxy ail d Pxx/Pzz were taken from two separated old 
measurement o 52 i 53 , such a comparison of the estimated 
experimental value of a l f- d with the calculated afj plus 
afj might not be very meaningful. Unfortunately, re- 
cent measurements of anomalous Hall and longitudinal 
resistivities on a same speciman were conducted only on 
polycrystalline Co metal o 55 ' 56 , and the deduced experi- 
mental values of the anomalous Hall conductivity a™ ,d 
are ~200 S/cm (Ref. 55) and -500 S/cm (Ref. M 
Clearly, further simultaneous measurements of both lon- 
gitudinal and Hall resistivities on the well-characterized 
thin films such as that on bee Fe reported in Ref. [l3| are 
needed to extract reliable anomalous Hall conductivities 
for hep Co, especially, a yz . 



B. Spin and orbital magnetic moments 

In addition to the AHC and SHC, including the SOI 
in the fully relativistic calculations would also give rise 
to nonzero orbital magnetic moments in the magnetic 



5 




r mk r A LH A 



FIG. 2. (color online) Relativistic band structures calculated 
with (solid lines) and without (dashed lines) on-site Coulomb 
interaction U = 1.6 eV. The magnetization is along the c-axis. 
The Fermi level is shifted to eV. 

solids concerned. The calculated spin and orbital mag- 
netic moments are listed in Table |TJ First of all, Table U 
shows that both the spin and orbital magnetic moments 
increases monotonically with increasing on-site Coulomb 
interaction U. However, this increase in the orbital mag- 
netic moment is much more significant than that of the 
spin magnetic moment. As mentioned earlier, due to its 
incomplete treatment of the 3d electron-electron correla- 
tion, the GGA calculations give an orbital moment which 
is only half of the corresponding experimental valued. 
In contrast, the theoretical orbital moment from the 
GGA+U calculations with U = 1.6 eV is almost doubled 
and becomes in rather good agreement with the experi- 
mental value (see the results for m//c in Table HJ. When 
the U is further increased, the orbital moment becomes 
much larger than the experimental value. This suggests 
that U — 1.6 eV would be appropriate for Co metals. 
Interestingly, as mentioned before, the orbital magnetic 
moments from the OPC calculations also agree well with 
the experiment o 22 ' 23 . 

Secondly, Table U indicates that the orbital magnetic 
moment depends significantly on the magnetization ori- 
entation, whereas the spin magnetic moment is not. Fur- 
thermore, this magnetocrystalline anisotropy in the or- 
bital magnetic moment increases pronouncedly when the 
on-site Coulomb interaction is taken into account. For 
example, the orbital moment anisotropy from the GGA 
calculations is about 5.0 % and becomes ~16.0 % in the 
GGA+U calculations with U = 1.6 eV. 



C. Effects of on-site Coulomb interaction on band 
structure 

To investigate how the on-site electron-electron cor- 
relation affects the electronic band structure and AHE 
in cobalt, we plot the relativistic energy bands along the 
high-symmetric lines in the Brillouin zone calculated with 



and without on-site Coulomb interaction for the magne- 
tization being along the c-axis in Fig. [2] and in the a-b 
plane in Fig. [3J The relativistic band structure of hep Co 
has been studied before by several researchers by using 
different band structure calculation methods (see, e.g., 
Refs. [H and [Til and references therein). The present 
relativistic GGA band structures (the red dashed curves 
in Figs. 2, 3 , and 4) are nearly identical to that reported 
in Refs. [H and 03- The relativistic band structure may 
be regarded as the results of a superposition of the cor- 
responding scalar-relativistic spin majority and minor- 
ity band structures with many accidental band crossings 
lifted by the SOI. These SOI-induced splittings in cobalt 
are, in general, much smaller than the exchange split- 
tings (see Fig. [5] and Fig [3]), and thus can be treated as 
a perturbation. 

Nonetheless, it is the inclusion of the SOI that makes 
the band structure depend on the magnetization orien- 
tation and this magnetization dependence of the band 
structure results in a strong anisotropy in such SOI- 
induced quantities as anomalous and spin Hall conductiv- 
ities. An example of clear changes in the band structure 
due to the rotation of the magnetization from the c-axis 
to an in-plane direction is the energy bands just below 
the Fermi level (Ep) in the vicinity of the H point (as 
highlighted by the dash-dotted green circle in Fig. 2 and 
Fig. 4a). These energy bands are significantly spin-orbit 
split when the magnetization is along the c-axis but they 
remain degenerate when the magnetization is in-plane 
(see Figs. 2, 3, and 4). Furthermore, the SOI-induced 
splittings for the bands near the Ep at the T point when 
the magnetization is in-plane, are nearly two times larger 
than that of the magnetization along the c-axis. These 
differences in the band structure near the Ep due to dif- 
ferent magnetization orientations thus give rise to a pro- 
nounced anisotropy in the SOI-induced phenomena such 
as the AHC, SHC and orbital magnetic moment. 

A pronounced change in the band structure due to the 
inclusion of the on-site Coulomb interaction is that all 
the d-dominant bands are lowered in energy by about 
0.2 eV relative to the Ep (see Figs. 2, 3, and 4). This is 
clearly seen in the energy window between -5.5 and -0.5 
eV. On the other hand, the energy bands near and above 
the Ep are much less affected when the on-site Coulomb 
interaction is taken into account, as shown in Figs. 2, 
3, and 4. In particular, for the in-plane magnetization, 
the energy bands in this energy range remain almost un- 
changed upon the inclusion of the on-site Coulomb in- 
teraction (see Fig. 3 and Fig. 4d). With respect to the 
AHC for the magnetization along the c-axis, the most 
significant change due to U is perhaps that the energy 
band in the vicinity of the T, L and H symmetry points 
(as highlighted by the green circles in Fig. 2 and Fig. 
4a) which is just below the Ep in the GGA calculations, 
is now pushed onto the Ep in the GGA+U calculations. 
As a result, the AHC <r£ y is increased by about 30 % 
upon including the on-site Coulomb interaction U = 1.6 
eV. 
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FIG. 3. (color online) Relativistic band structures calculated 
with (solid lines) and without (dashed lines) on-site Coulomb 
interaction U = 1.6 eV. The magnetization is in the a-b plane. 
The Fermi level is shifted to eV. 

In Fig. IU we display the calculated anomalous Hall 
conductivity (c^) and the number of valence electrons 
per Co atom(n e /atom) as a function of the Ep, together 
with the relativistic band structures without and with on- 
site Coulomb interaction using U — 1.6 eV and J = 0.9 
eV. The fine k point mesh of 130x130x69 is used. Fig. 
3b shows that in general, the a~y calculated with and 
without the on-site Coulomb interaction are rather sim- 
ilar. The most significant difference is that the a^ y at 
the true Fermi level (0 eV) from the GGA+U calcula- 
tion is 30 % larger than that obtained from the GGA 
calculation. As mentioned in the preceding paragraph, 
this enhancement of the cr^y could be attributed to the 
fact that the inclusion of the on-site Coulomb interaction 
pushes the energy band near the T, L and H points to 
just above the Fermi level (Fig. 2 and Fig. 4a). Fig. 

indicates that for m//ab, the two cjy Z spectra calcu- 
lated with and without the on-site Coulomb repulsion are 
nearly identical, except their magnitudes differ slightly 
here and there. 

Fig. 4b also shows that the calculated a^ y are quite 
flat in the vicinity of the true Fermi level. As the Fermi 
level is artificially lowered from about -0.3 eV to -0.71 eV 
(n e /atom = ^7.95), the a^ y increases steadily and even- 
tually peaks at -0.71 eV with a maximum value of ^2600 
S/cm. When the Fermi level is further lowered, the a^y 
first oscilates and then starts to decrease sharply at -1.2 
eV, and eventually, it changes its sign at -1.3 eV. If the 
Fermi level is artificially raised above eV, the a^y de- 
creases gradually and then changes to the negative value 
at ^0.5 eV. The reaches a maximum magnitude (- 
1900 S/cm) at ~1.15 eV. Above this energy, the mag- 
nitude of the a£ decreases steadily and becomes very 
small above ~ 2.0 eV. Fig. 4e suggests that the calcu- 
lated <jy Z for raj jab is overall similar to the a£ for m//c 
(Fig. 4b), especially above -1.5 eV. In particular, the a yz 
spectra also have a negative peak (-800 S/cm) near 1.25 
eV and a positive peak (1500 S/cm) around -0.70 eV. 




a A (S/cm) n e (e/atom) 



FIG. 4. (color online) Relativistic band structures calculated 
with (solid lines) and without (dashed lines) on-site Coulomb 
interaction U = 1.6 eV (a, d), anomalous Hall conductivity 
(<7ij) (b, e), and number (n e ) of valence electrons per Co atom 
(c, f). Upper panels are for the magnetization along the c-axis 
(to/ /c) and lower panels are for the magnetization in the a-b 
plane (to/ /ab). The Fermi level is shifted to eV. The green 
circles in (a) are a guide to the eyes only. 

Nonetheless, the peaks in the a yz are smaller than the 
corresponding ones in the a^ y spectra. 

D. Spin polarization of Hall current 

Electric and spin current transports are determined 
by the characteristics of the band structure near the 
Fermi level in the solids concerned. For the spintronic 
applications, it would be interesting to examine the ex- 
tent that the density of states (DOS) near the Ep of 
the systems considered are spin-polarized. This spin- 
polarization (p D °S) of the DOSs near the Ep is defined 
as 

pDOS = N t (E F ) - N^Ep) 

N^^ + N^Ef) 1 1 1 

where Nf(Ep) and N±(Ep) are the spin-up and spin- 
down DOSs at the Ep, respectively. p DOS mav vary 
from -1.0 to 1.0. For the half-metallic materials, P DOS 
is either 1.0 or -1.0, and the charge current would be 
fully spin polarized. We have calculated DOSs of hep 
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Co with the magnetization being along the c-axis and in 
the a-b plane, and the calculated spin-decomposed DOSs 
at the Ep are listed in Table ITT1 The corresponding spin 
polarization p DOS defined above for both magnetization 
directions is also listed in Table [TTJ It is clear from Table 
HI] that the spin polarization is negative and the degree of 
spin polarization is rather high, being 66 ~ 69 %. Fur- 
thermore, the spin polarization p DOS is independent of 
magnetization orientation, in strong contrast to the high 
anisotropy of the calculated anomalous and spin Hall con- 
ductivities (Table J}. We note that the current spin polar- 
ization of cobalt measured in the tunneling experiments 58 
and in the Andreev reflection experiments^ 9 - is about 42 
% and 45 %, respectively, being considerably lower than 
the calculated p DOS . 

As pointed out by several researchers before (see, e.g., 
Ref. [6fJ), the spin polarization p DOS defined by Eq. 
(3) is not the spin-polarization of the transport cur- 
rents measured in the experiments. The anomalous Hall 
effect has recently received intensive renewed interest 
mainly because of its close connection with spin transport 
phenomenal. In particular, for spintronic applications, it 
could be advantageous to use the Hall current from ferro- 
magnets as a spin-polarized current source, instead of the 
longitudinal current, because of the topological nature of 
the intrinsic AHE 27 . Therefore, it would be interesting to 
know the spin polarization of the Hall current. To this 
end, let us define the spin-polarization Pff of the Hall 
current as 



pH _ a ij 



Hi 



<7, 



Hi > 



(4) 



where and ov- are the spin-up and spin-down Hall 
conductivities, respectively. Note that this decomposi- 
tion of the Hall conductivity in terms of the simple two 
current model, works reasonably well only for metals con- 
taining light elements such as 3d transition metals, but 
may fail for metals containing heavy elements such as 
Coi-zPt^ alloysSi. The af^ and af^ can then be ob- 
tained from the calculated AHC a{\ and SHC af, via the 
following relations 



iff 



HI 



HI 



(5) 



The 0"jj , o^j- , and P?j obtained in this way, are listed 
in Table |Hj Note that, unlike the spin-polarization of 
the DOSs and also of the longitudinal charge currents, 
the magnitude of the P?j can be larger than 1.0. And 
this is because spin-decomposed Hall currents can go ei- 
ther right (positive) or left (negative). In the nonmag- 
netic materials, the charge Hall current is zero, and hence 
&ij = — ffy . In this case, we would then have pure spin 
current, and PM = oo. 

Interestingly, TablelHlshows that, in contrast to PP , 
the spin polarization of charge Hall current is highly 
anisotropic. In particular, when the magnetization is 



parallel to the c-axis, both spin-up and spin-down Hall 
conductivities are positive and hence the magnitude of 
Pfj is smaller than 1.0. The magnitude of both P^ for 
U = and 1.6 eV is rather large, being 0.48 and 0.61, re- 
spectively, although it is smaller than the corresponding 
pDOS value. In contrast, for the in- plane magnetiza- 
tion with [7 = and 1.6 eV, the spin- up and spin-down 
Hall conductivities have an opposite sign. Therefore, the 
magnitude of P& is larger than 1.0, being 4.90 and 1.59 
for Pf- for U = and 1.6 eV, respectively. Remark- 
ably, this suggests that the Hall current could be a good 
spin-polarized current source. The best materials known 
for the spin-polarized current source are pressumably the 
half-metallic magnets with the spin-polarization being 
1.0 (see, e.g., Ref. [62| and references therein). The lon- 
gitudinal spin-polarized current from a magnetic metal 
is charge current while the Hall current could consist of 
both spin-polarized charge current and pure spin current, 
as mentioned before. Pure spin current is dissipationless 
while charge current would consume energy— Therefore, 
using the Hall current instead of the longitudinal cur- 
rent as the spin-polarized current source could also save 
energy. 

Table iHl also indicates that Pf- can be significantly af- 
fected by the on-site Coulomb interaction. For example, 
upon the inclusion of the on-site Coulomb interaction 
U = 1.6 eV, the magnitude of Pff would be increased 
by 25 % when the magnetization is along the c-axis, but 
would be decreased by ~ 65 % when the magnetization 
is in the a-b plane (see Table HH) . There has been no re- 
port on the measurement of the spin polarization of the 
anomalous Hall current. Therefore, one could be tempted 
to compare the theoretical P& with the measured spin 
polarization (P L ) of the longitudinal charge curren t 58 ' 59 . 
However, the theoretical P^ z calculated using U = and 
1.6 eV is very different from the measured P L , although 
the theoretical P^ v is comparable to the P L (see Table 
[TTJ) . This suggests that the spin polarizations of the lon- 
gitudinal and Hall currents are two different things and 
may have no correlation. It is hoped that the interesting 
findings reported in this paper will stimulate the mea- 
surements on the spin polarization of the Hall current in 
either hep Co or other magnetic metals. 



IV. CONCLUSIONS 

We have presented theoretical intrinsic anomalous and 
spin Hall conductivities, spin and orbital magnetic mo- 
ments and also spin polarization of Hall currents in hep 
Co with different magnetization orientations calculated 
in both the GGA and GGA+J7 schemes. The accu- 
rate FLAPW method is used. First of all, we find that 
the AHC (v^y) and orbital magnetic moment (m ) of 
cobalt with m//c calculated in the GGA+U scheme with 
U = 1.6 eV and J = 0.9 eV are in good agreement with 
the corresponding experimental values while, in contrast, 
the cr^y and m from the GGA calculations are signifi- 
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TABLE II. Theoretical spin-decomposed Hall conductivities 
[(cr^ and (S/cm)] and spin-polarization of the Hall cur- 
rent (P H ) as well as spin-decomposed densities of states (Nf 
and N_i) (states/eV/unit cell) and spin-polarization of the 
density of states at the Fermi level (P DOS ). The on-site 
Coulomb interaction U = 1.6, 1.9 and 2.5 eV are used, with J 
= 0.9 eV being constant. The measured longitudinal current 
spin-polarizations (P L ) are also listed for comparison. 

GGA GGA+U GGA+U GGA+U Expt. 
17 = 1.6 eV U = 1.9 eV U = 2.5 eV 









m| c 




^ xy 


125 


121 


127 


96 


a Al 

w xy 


359 


497 


517 


732 


pH 

1 xy 


-0.48 


-0.61 


-0.61 


-0.77 


N t 


0.316 


0.281 


0.291 


0.257 


Nr 


1.485 


1.509 


1.504 


1.400 


pDOS 


-0.66 


-0.68 


-0.68 


-0.69 


pL 











-0.45 b 









m| ab 






-11.9 


-11.5 


-6.0 


11.3 


a Al 
u yz 


18.0 


50.4 


59.8 


70.5 


pH 

± yz 


-4.90 


-1.59 


-1.22 


-0.72 


Nf 


0.302 


0.293 


0.290 


0.253 




1.491 


1.539 


1.554 


1.626 


pDOS 


-0.66 


-0.68 


-0.68 


-0.69 



a Tunneling experiments (Ref. \E3tl . 

b Andreev reflection measurements (Ref. [591) , 



cantly smaller than the measured ones. This suggests 
that the on-site Coulomb interaction, though moderate 
only in Co metals, should be taken into account properly 
in order to get SOI-induced properties. The most obvi- 



ous effect of including the on-site Coulomb interaction on 
the band structure is that all the d-dominant bands in 
the energy window between 0.5 and 5.5 eV below the Ep 
are lowered in energy relative to the Ep by about 0.2 eV 
while the energy bands above the Ep remain nearly unaf- 
fected. In the context of the Hall conductivities, the most 
significant change due to U is that the energy band in the 
vicinity of the T, L and H symmetry points just below the 
Ep in the GGA calculations, is pushed onto the Ep in 
the GGA+{7 calculations, resulting in a ~30 % increase 
in the a^ y and also a ~60 % increase in the af y . Sec- 
ondly, we find that the calculated AHC and SHC (cr s ) are 
highly anisotropic, and the ratio of the Hall conductivity 
for m//c to that for m/ /ab from the GGA+C/ calcula- 
tions, is -16.0 for the AHC and -6.0 for the SHC. Fi- 
nally, the calculated Hall current spin polarization (P H ) 
is -0.61 for to/ /c and -1.59 for m/ /ab. The theoretical 
a s and P H are also found to be significantly affected 
by the inclusion of the on-site Coulomb interaction U. 
We hope that this work would stimulate simultaneous 
measurements of both longitudinal and Hall resistivities 
on the well-characterized thin films to extract reliable 
experimental anomalous Hall conductivities for hep Co, 
especially, <jy Z . 
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